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Transversal interface dynamics of a front connecting a stripe pat- 
tern to a uniform state 
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Abstract. - Interfaces in two-dimensional systems exhibit unexpected complex dynamical be- 
haviors, the dynamics of a border connecting a stripe pattern and a uniform state is studied. 
Numerical simulations of a prototype isotropic model, the subcritical Swift-Hohenberg equation, 
show that this interface has transversal spatial periodic structures, zigzag dynamics and com- 
plex coarsening process. Close to a spatial bifurcation, an amended amplitude equation and a 
one-dimensional interface model allow us to characterize the dynamics exhibited by this interface. 



Introduction. — Non equilibrium processes often 
lead in nature to pattern formation developing from a ho- 
mogeneous state through the spontaneous breaking of the 
symmetries present in the system [1]. In recent decades, 
much effort has been devoted to the study of pattern for- 
mation (see review [2] and the references therein) arising in 
systems such as chemical reactions, gas discharge systems, 
C0 2 lasers, liquid crystals, hydro dynamic or electrocon- 
vective instabilities, and granular matter (see review [3]), 
to mention a few. A unified description for the dynam- 
ics of spatially periodic structures, developed at the on- 
set of a bifurcation, is achieved by means of amplitude 
equations for the critical modes. Such a description is 
valid in the case of weak nonlinearities and for a slow spa- 
tial and temporal modulation of the base pattern [2], As 
an example, the Newell- Whitehead-Segel equation [4] de- 
scribes the dynamics of a stripe pattern formed in a two- 
dimensional system. Another ubiquitous phenomenon in 
nature is the interface dynamics or front propagation. The 
concept of front propagation, emerged in the field of pop- 
ulation dynamics [5-7] , has gained growing interest in bi- 
ology, chemistry, physics, and mathematics (See, e.g. [8] 
and references therein). These interfaces connect two ex- 
tended states, such as: uniforms states, patterns, oscilla- 
tory, standing waves, spatio-temporal chaotic and so forth. 

In one-dimensional systems, an interface connecting two 
uniform stable states, the most favorable state — for in- 
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stance energetically — invades the other one with a con- 
stant and unique speed [9,10]. This speed is zero, that 
is the front is motionless, at the Maxwell point [11]. The 
above picture changes, when one considers an interface 
connecting a pattern state and a uniform one or two pat- 
terns. Due to spatial translation symmetry breaking, the 
interface is motionless in a range of parameters, the pin- 
ning range [11-13]. This behavior is called locking phe- 
nomenon or pinning effect. In bidimcnsional dynamical 
systems, few experimental and theoretical studies have 
been performed on fronts connecting patterns and uniform 
states [14-17]. 

The aim of this letter is to study the dynamical behav- 
iors of a front connecting a stripe pattern to a uniform 
state. Numerical simulations of a prototype model — the 
isotropic Swift-Hohenberg equation — show that the lock- 
ing phenomenon of a flat interface persists. However, this 
flat interface is nonlinearly transversely unstable, that is, a 
finite perturbation of this interface leads to the appearance 
of an initial wave number which is subsequently replaced 
by zigzag dynamics, which presents a complex coarsen- 
ing. Increasing the longitudinal interface size, the flat 
interface exhibits a transversal spatial instability, which 
originates a periodical structure at the interface. We have 
termed these interfaces embroideries. In order to explain 
these behaviors in an unified manner, we make use of the 
amended Newell- Whitehead-Segel equation and prototype 
one-dimensional model for the interface, which describe 
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properly the dynamics observed on the interface connect- 
ing a stripe pattern with a uniform state. 

The model. — A simple isotropic model which ex- 
hibits coexistence between a uniform state and a stripe 
pattern is (the subcritical Swift-Hohenberg equation [2]) 

d t u = eu + vu z - u 5 - (V 2 + q 2 ^) 2 u, (1) 

where u (x, t) is an order parameter, e is the bifurcation 
parameter, q is the wave number of the stripe pattern, 
v the control parameter of the type of bifurcation (su- 
percritical or subcritical), and V 2 is the Laplacian oper- 
ator. The model (1) describes the confluence of a sta- 
tionary and a spatial subcritical bifurcations with reflec- 
tion symmetry (u — > —it), when the parameters scale as 
u ~ £ 1/4 , v ~ e 1 / 2 , q ~ e 1 /*, d t ~ e and V ~ e 1 / 4 
(e -C 1). The above model is often employed in the de- 
scription of patterns observed in Rayleigh-Benard convec- 
tion [2]. For small and negative e and — 9^ 2 /40 = e sn < 
e < 0, the system exhibits coexistence between a uni- 
form state u(x, y,i) = and a stripe pattern u (x, y, i) = 

y/V (^2(1 + yi + 40e/9z/) cos (<f • f)^ + o(^ 5 / 2 ), where q 

is an arbitrary vector with modulus q. For e = e sn the 
model has a saddle-node bifurcation that give rises to sta- 
ble and unstable patterns. For e — the uniform state 
becomes unstable. When one considers the above model 
in one spatial dimension, it is well-known that in the co- 
existence region (e sn < e < 0) the model exhibits a front 
connecting a spatially periodic solution and uniform state. 
For e close to (e sn ), the pattern (uniform) state in- 
vades the uniform (pattern) state, i.e. the system displays 
front propagation. The front is motionless at the inter- 
val S- < e < £+, the pining range [11]. Note that, in 
this range, the state which is energetically more favorable 
does not invade the less favorable one. It is important to 
note that the interface has two characteristic lengths, the 
pattern length 2ir/q and the interface size, which is rep- 
resented by A in the inset of Fig. ??. Decreasing v the 
interface size increases. 

Numerical results. — In two spatial dimension the 
above scenario changes drastically. If one considers a sim- 
ilar parameters setup (e sn < £_ < e < e + < 0), the flat 
interface of model (1) exhibits locking phenomenon (cf. 
Fig. ??a). However, when the longitudinal interface size A 
is increased the interface suffers a supercritical transversal 
spatial instability and it gives rise to transversal periodic 
structures. The typical observed structures are depicted 
in Fig. ??. We term these periodical structures at the in- 
terface embroidery. The embroidery size is proportional to 
the longitudinal interface size. The embroideries are con- 
sequence of two facts: the spatial isotropy and the inter- 
action of enveloped variation with the underlying pattern 
(pinning effect). Seeing that at the interface the stripe 
pattern can develop in any direction as consequence of 



the isotropy, however the interaction of enveloped vari- 
ation with the underlying pattern freezes the growth of 
these unstable modes, as we shall see with the amended 
amplitude equation. Based on the ID theory of front in- 
teraction which predicts a family of localized patterns [18] , 
one expects to find a family of localized stripes, that is, the 
localized stripes are the transversal expansion of ID local- 
ized patterns. The typical stable localized states observed 
in model (1) are illustrated in Fig. 2. It is important 
to note that, localized stripe patterns with flat interfaces 
have been observed in ref. [19] . 

Numerically, we observe that in a finite region of param- 
eters the embroidery interface is linearly stable, i.e. this 
interface with a small initial perturbation evolves to the 
interface without perturbations. However, when we con- 
sider a large perturbation the interface does not evolve to 
the interface without perturbation. This interface exhibits 
a nonlinear zig-zag instability, which is characterized by a 
transversal instability without a well-defined wavelength. 
In spite of this, initially this instability has a well defined 
wave number close to 2n/q, later on, the sinusoidal inter- 
face becomes an angled line composed of pieces of interface 
turned with well define angles, zig- facet or zag- facet. It is 
important to note that the flat interface-in the regime of 
parameters where it is linearly stable — exhibits the same 
behavior, that is, a finite initial perturbation of the flat 
interface gives rise to a zig-zag dynamics. In order to illus- 
trate this nonlinear instability, Fig. ??a shows two embroi- 
dery interfaces, one with a initial small perturbation (left 
interface) and the other with a finite one. After a finite 
time, the interface with initial small perturbation evolves 
to an embroidery interface. While the Right-interface de- 
velops a zig-zag structure and the pattern propagates over 
the uniform state (cf. Fig. ??b). Hence, the flat and em- 
broidery interfaces are nonlincarly transversely unstable. 

In the zig-zag interface, two adjacent facets whose ori- 
entations are opposite, are connected by a region of strong 
curvature that we term a corner. The dynamics shown by 
the zig-zag interface consists then in reassembling domains 
of even orientation, the angle facets staying unchanged, 
which is a coarsening dynamics. This process occurs due 
to annihilations of corners and without a characteristic 
length scale. Actually, the averaged domain size increases 
regularly in time. Simultaneously to this coarsening pro- 
cess, the zigzag interface propagates from one extended 
state to the other one. Figure 4 depicts the typical coars- 
ening process observed at the interfaces. 

The amended amplitude equation. To explain 
the dynamics exhibited by the interfaces of model (1) in a 
unified framework, close to the spatial bifurcation (e <C 1, 
and v ~ e 1 / 2 ), we can introduce the ansatz 

u(x, y,t) =a A(x = f , y = -E=, r = pi e** 
8q 2 24q 2 
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where a = l = 2qVW/3y/\s\, t = 

10/9^ 2 |e|, and A(X,Y,t) is the envelope of the pattern 
that describes the front solution (when the envelope is 
uniform and not null the initial model (1) has an stripe 
pattern in the ^-direction) , the "• • •" stands for the high 
order terms in the envelope A, c.c. means complex conju- 
gate, and {X,Y,t} are slow variables. In this ansatz (2), 
we consider that q is of order one or larger than the other 
parameters. Introducing the above ansatz in equation (1), 
we find the following solvability condition for the envelope 
(Amended subcritical Newell- Whitehead-Segel equation) 

d T A = je + \A\ 2 \Af + [d x ^) 2 1 A + r,A 3 e™, 

(3) 

where e = -lOe/9^ 2 , k = 2Vl0q/3uy/e, and r, = 1/3 for 
model (1). Using the above scaling the terms inside brack- 
ets are order e [20] , and the last term is exponentially small 
which is regularly neglected in the multiple scaling ap- 
proach. However, to account for the coupling between the 
large scale envelope and the small scale underlaying the 
stripe state, we consider n as free parameter and k a finite 
number. Notice that consider other solvability conditions 
for the envelope one can obtain diverse small exponential 
terms for equation (3) of the form A n A m e lq ^ n ~ 1 ~ m ^ x . For 
the sake of simplicity we have consider the dominant one in 
A, which correspond to the last term in Eq. (3). Neverthe- 
less, we expect that these entire exponential small terms 
have qualitatively the same effect in the dynamics. When 
the term proportional to spatial forcing is zero (77 — *- 0), 
the above model is the Newell- Whitehead-Segel amplitude 
equation, which have been used deeply to explain the ap- 
pearance of stripe pattern [2]. Nevertheless, this model 
does not account for locking phenomenon [11]. The inclu- 
sion of spatial forcing terms in the amplitude equations 
in ID allows understanding the locking phenomenon, the 
pinning range and localized structures [12,13,18]. In the 
extreme limit, 77 — > 0, it is straightforward to show that 
model (3) has a front solution connecting two homoge- 
neous states, and (l + \J\ + 4e) /2, when e < 0. Which 
accounts for the connection between the stripe pattern 
and the uniform state. These two homogeneous states 
are energetically equivalent at eu — —3/16 — the Maxwell 
point — and it has the form 

A=J 3 /i e« (4) 

where P stands for the position of flat interface, and 9 
is an arbitrary phase. This flat interface is transversally 
unstable and it gives rise to a complex coarsening pro- 
cess, zigzag instability. Numerically we have computed 
the growth rate of each mode with wave number k of the 
flat interface — spectrum. Fig. ??a depicts this instability 
and the spectrum of the flat interface. Note that although 
the numerical simulations have been done at the Maxwell 



point, the interface propagates from a state that repre- 
sents the stripe to a uniform one. When 77 is small the 
non-null uniform solution becomes a stripe state in the y- 
direction (cf. Fig. ??b). The zigzag dynamics exhibited 
by the interface disappears and is replaced by a transver- 
sal periodic structure, embroidery (cf. Fig. ??b). Hence, 
these transversal embroidery structures are consequence 
of the interaction of the spatial forcing — generated by un- 
derlying pattern — with the transversal instability of the 
Newell- Whitehead-Segel model. Changing e the embroi- 
dery interface is motionless in a range of parameter of e, 
the pinning range. Fig. ??b shows the typical motionless 
embroidery interface observed in model (3) and the re- 
spective spectrum of the flat interface. Increasing 77 the 
embroidery interface becomes a flat interface and the am- 
plitude of the stripes increase. Fig. ??c shows the stable 
flat interface and its respective spectrum. The flat and 
embroidery interfaces exhibit by the amended amplitude 
equation are nonlinear unstable, that is, a finite perturba- 
tion of these interfaces give rise a zigzag dynamics. There- 
fore, the dynamics exhibits by the universal model (3) is 
similar to those shown by prototype model (1). 

Phenomenological ID model. The model (3) 
seems simpler than Eq. (1), however, the analytical de- 
scription of the interfaces in these model is a thorny task. 
A standard method to grasp the dynamics exhibited by 
the interface of the precedent models is to derive a one- 
dimensional equation for it. This method consists in using 
as an ansatz, the front solution (5) plus a small correction, 
that is, 

A=\J ft + w Q (X-P,Y)\e i9 , 

(5) 

where the continuous parameter P is promoted to a field 
(P(Y,t)), wq is a small complex correction function, 
which is order of spatial variation of the position of the 
interface (wq <~ dyyP) [21-24]. However, we can not use 
this weakly nonlinear method because the more unsta- 
ble transversal mode has a finite wave number (q). This 
type of method requires that unstable modes have small 
wave number. To understand the mechanism of the dif- 
ferent structures observed at the interface, and based on 
symmetry arguments for the interface [21-24] and on the 
effect of spatial forcing term [12, 13], we propose the fol- 
lowing phenomenological equation for the position of the 
interface (convective and forced Cahn-Hilliard equation) 

d T P = ePyy + P y Pyy - Pyyyy + aPy - k sin(QP) + A 

(6) 

this model has the equilibria P n = mr/Q, n = 0, ±1, 
±2, • • •. The equilibria Pim as function of forcing k have 
similar spectrum to those shown in Fig. ??. For large k, 
these equilibria are stable, i.e. this model has a family 
of stable flat interfaces. Decreasing n these flat interfaces 



p-3 



Marcel G. Clerc 1 Daniel Escaff 2 Rene Rojas 2 1 



become unstable and give rise to spatial periodical state, 
embroideries. Finally, for small {a, k} this model exhibit 
zigzag instability characterized by logarithmic and power 
law for the coarsening. Hence, the one-dimensional model 
(6) presents qualitative similar dynamics that those shown 
by models (1) and (3). In the coexistence region of uniform 
states, one expects to find stable localized horm solution 
[7], and then we expect to find localized horm state at 
the interface. Fig. ?? illustrates the typical stable horm 
solutions of one-dimensional model (6) and the prototype 
model (1). The study of the effect of the forcing on zigzag 
dynamics of above model is in progress. 

Summary. — Isotropic systems which have coexis- 
tence between stable stripe pattern and uniform states 
can exhibit interfaces connecting these states. These in- 
terfaces present a rich and unexpect transversal dynam- 
ics like transversal spatial instability, nonlinear zig-zag in- 
stability and localized states. Recently in Ref. [25], it is 
shown that in anisotropic systems the flat interface link- 
ing rolls pattern with uniform one is transversal stable. 
Hence, the wealthy transversal dynamics is a consequence 
of the spatial isotropy and the pinning effect. 

Although the prototype model under study Eq. (1) is 
variational — the dynamical evolution of this model has the 
tendency to minimize its Lyapunov functional — the inter- 
face dynamic exhibited by this model is robust. Since, 
this interface dynamics is well described by the amended 
amplitude equation, which is valid close to the spatial in- 
stability. Hence, we expect that a system that exhibits a 
coexistence between a stripe pattern and a homogeneous 
state should be present a rich interface dynamics. 
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